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0.   ABSTRACT 

Structures  are  subject  to  changing  loads  from  various  sources. 
In  many  instances  these  loads  fluctuate  in  time  in  an  apparently 
random  fashion.   Models  are  considered  for  which  the  stress  put 
on  the  structure  by  various  loads  simultaneously  can  be  described 
by  a  regenerative  process.   The  distribution  of  the  first  time 
until  the  stress  on  the  structure  exceeds  a  given  level  x  is  studied 
Asymptotic  properties  of  the  distribution  are  given  for  a  large 
stress  level  x  and  for  the  tail  of  the  distribution  for  fixed 
finite  stress  level  x.   Simulation  results  are  given  to  assess 
the  accuracy  of  using  the  asymptotic  results  to  approximate  the 
distribution . 


1.   INTRODUCTION 

Many  physical  structures  are  subject  to  varying  physical 
loads  from  various  sources:   wind,  snow,  and  earthquakes  being 
examples.   In  many  instances  the  total  load  experienced  by  a 
structure  varies  in  time  in  an  apparently  random  fashion.   Certain 
load  components,  e.g.  those  of  snow  and  ice  accumulation,  vary 
rather  slowly;  others  such  as  those  associated  with  winds  or 
earthquakes  occur  more  nearly  as  impulses.   The  problem  is  to 
design  structures  to  withstand  a  coincidence  of  such  loads  with 
approximately  a  prescribed  (high)  probability.   In  engineering 
terms  we  wish  to  work  towards  developing  a  rational  safety  factor 
criterion  for  designing  structures  to  withstand  the  combination 
of  loads  anticipated. 

In  this  paper  we  will  study  the  time  for  the  load  combination 
process,  (X(t);  t  >_  0 }  to  exceed  a  given  stress  level  for  the  case 
in  which  the  load  combination  process  is  a  regenerative  process 
taking  non-negative  values  (c.f.  Cinlar  [1975],  page  298).   Many 
of  the  models  for  load  combinations  that  have  been  studied  are 
regenerative  processes;  (cf.  Pearce  and  Wen  [1983],  Wen  [1981], 
Shanthikumar  and  Sumita  [19  83]) . 

EX/AMPLE  1.   The  load  combination  process  is  the  superposition 
of  two  load  types:   shock  loads  and  constant  loads.   A  constant 
load  (e.g.  caused  by  snow  or  rain  accumulation)  exists  at  one 
level  for  a  random  time  and  then  changes  to  a  new  level.   Let 
Y(t)  denote  the  constant  load  magnitude  at  time  t.   The  instants 
of  change  in  the  magnitude  of  Y(t)  occur  according  to 


a  renewal  process  with  inter-renewal  distribution  H;  successive 
magnitudes  of  the  constant  load  process  are  independent  iden- 
tically distributed  with  distribution  F.   Impulse  or  "shock" 
loads  (e.g.  caused  by  wind  gusts  or  earthquakes)  occur  at 
random  moments  and  have  time  duration  of  length  0.   Let  Z(t) 
denote  the  shock  load  magnitude  at  time  t.   Given  the  constant 
load  magnitude  at  time  t,  Y(t)  =  y,  the  probability  a  shock 
load  will  occur  in  the  time  interval  [t,t+h]  is  y(y)h  +  o(h); 
the  magnitude  of  the  shock  load  is  conditionally  independent 
of  everything  given  Y(t)  =  y  and  has  distribution  G(y;*) . 

The  load  combination  process  magnitude  at  time  t  is 
X(t)  =  Y(t)  +  Z(t).   {X(t)  ;  t  >_  0}  is  a  regenerative  process 
with  regeneration  times  the  times  of  change  of  the  constant 
load  process. 

Gaver  and  Jacobs  [1981]  studied  a  special  case  of  the 
above  model  in  which  the  interrenewal  distribution  H  is  exponen- 
tial, and  the  shock  load  process  is  a  compound  Poisson  process 
independent  of  the  constant  load  process  with  shock  arrival 
rate  y  and  magnitude  distribution  function  G.   Other  similar 
load  combination  models  have  been  studied  in  the  past;  cf. 
Peir  and  Cornell  [1973] ,  Wen  [1977] ,  Pearce  and  Wen  [1983] . 

Let  X(t)  denote  the  magnitude  of  the  load  combination 
process  at  time  t.   The  process  {X(t)  ,  t  >_  0  }   is  assumed  to 
be  a  regenerative  process.   Let 


T   =  inf {t  >  0:  X(t)    x}  ,  (1.1) 


the  first  time  the  load  combination  exceeds  a  level 


x 


Section  2  is  concerned  with  asymptotic  distribution  of  T 
as  x  ■*■   oo.   It  is  shown  that  under  certain  assumptions  the  distri- 
bution of  the  normalized  random  variable  T  (E[T  ]~  )  is  approximately 
unit  exponential  for  large  x,  and  error  bounds  on  the  rate  of 
convergence  are  obtained.   This  result  is  related  to  that  of 
Keilson  [1979,  page  134]. 

In  Section  3  the  tail  of  the  distribution  of  T   for  finite  x, 

x 

P{T  >t},  will  be  studied.   An  asymptotic  result  concerning  the 
exponentiality  of  P{T  >t}  for  large  t  is  given. 

In  Section  4  simulation  results  are  presented  to  study  the 
accuracy  of  using  the  two  .  asymptotic  results  to  approximate  the 
distribution  of  Tx . 


2.   THE  DISTRIBUTION  OF  THE  FIRST  PASSAGE  TIME 
FOR  THE  LOAD  COMBINATION  PROCESS 

Let  {X(t);  t>^0}  be  a  regenerative  process  taking  non-negative 

values  representing  the  load  combination  process;  [cf .  Cinlar  [1975]  , 

page  298] .   Let  S   denote  the  nth  regeneration  time;  {S  }  is  a 
^  ^  n  ^  n 

renewal  process.   We  will  assume  there  is  a  regeneration  at  time 

0.   Let  T   be  as  in  (1.1)  and  put  c|>(£)  =  E[exp(i^T  )],  the  Fourier 
x  x 

transform  of  T  .   A  renewal  theoretic  argument  yields 
x 


(J)(C)   =   E[exp(UTx};Tx  <  S-^  +  E  [exp  { i?Tx }  ;  S±    <    T^\  (2.1 

=   E[exp{i£T  };T     S,  ]  +  E  [exp{  i£S,  } ;  S,  <  T  ]<(>(£) 

Thus, 


i£T 

E[e         X;T      <    S,  ] 

♦  (5)    =  ^sf •  (2'2> 

1  -E[e         1;S1    <    Tx] 


Similarly,    if   m(x)    =   E[T    ],    then 


m(x)     =    E[T    ;T  S, ]     +    E[T    ;    T       >    Sn  ]  (2.3) 

X  X  X  XXX 


E[T    ;T       <    S, ]     +    E[S, ;T       >    S, 
x      x   —      1  1      x  1 


+    P{T      >    S, }m(x) 

X  X 


Therefore, 


E[min(Tx,S1)  ]  E  [min  (Tx  ,  S^  ] 

mU)     =    1  -P(T      >    S,}       =  PIT      <    Sl>  '  (2'4) 

x  1  x  —      1 


The    following   assumptions   will    be   made    for    the   remainder 
of    this    paper. 


0    <    E[S2]     <    -       ;  (2.5) 


lim  P{T      <    St }    =    0     .  (2.6) 

x   —      1 

X->oo 


It   now    follows    from    (2. 4) -(2. 6)     that 


lim  m(x)    =   °°    .  (2.7) 

x-*-°° 


LEMMA    1 


TTI"1*'^-  r^icl  E[si;Si  <  V 

E   1  2  1/2     1/2 

i.    g(X)  2E[min(Sl,Tx)]  +  E[s2]i/26(x)1/2 

+  UlB(x)3/2     E[Sl'1/2 


2E[min(S, ,T  ) ] 
1   x 


where  3 (x)  =  P{T   <  S,  } . 

x  —   1 

Proof.   It  follows  from  (2.2)  that  there  are  random  variables 
e,  and  £„  such  that 

4>(— £v)  =  E[l  +  -4^T  T   +  e_;T   <  S,  ]  (2.8) 

y  m(x)  m(x)   x     1   x  -   1 


x    [1    -    E[l    +   5f|rS1    +    e2;Sl    <    TJ]"1 


[P{Tx   1  Sl}    +   mTxT  E[VTx   1  Sl]     +    E[el'Tx   1  Sl 


x  [p{Tx  i  si}  "  mrir  E[si;Si  <  Tx]  "  E[£2;Si     Tx]]_1 


=    (m(x)P{T      <    S,  }  +  if,E(T    ;T       <    S,  )    +   m(x) E [ e, ; T  S,  ]  ) 

x   —      1  xx—      1  1      x    —      1 


[m(x)P{T       <    S1}  -i?E[S1;S1    <    Tx]     -    m(x)E[£2;S1    >    Tx]  ] 


(a   +    i£b)    x  (c   -    i£d)       . 


The    difference    between    the   Fourier    transform   of   T    (E[T    ] 
and    that    for   a   unit   exponential    is 


[(a  +i£b)    x  (c  -i£d)    4-    (1  -  i£)    1  (2.9) 


=    [(a-c+C2b)    +   i£(b-a+d)]    x  [ (c-i£d) (l-i£) ]    1 


_    NUM 
"    DEN 


The    term 


a-c+C2b   =  m(x)P{T      <    S,  }    +  m(x)  E  [e,  ;T      <_  S,  ]  (2.10) 

-    m(x)P{T      <    Sn }    +   m(x)E(£n;Sn     <    T    ] 
x    —      1  2       1  x 

+    C2E[T     ;T       <    S, ] 
x      x   —      1 


m(x)([E[£;L;Tx    <    S±]     +    Ef^,^    <    Tx]  ) 


+    ?2E[Tx;Tx    <    S1] . 


It    follows    from    (2.4)     that    the    term 


b-a+d   =    E[T    ;T      <    Sn  ]     -     [m(x)P{T       <    S,}    +    m(x)E[e,;T       <    Sj] 
xx—      1  x—      1  1      x   —      1 


+    E[S     ;S,     <    T    ]  (2.11) 

J_  J-  X 


b-a+d    =    E[min(T    ,S    )]     -    E[min(T    ,S,)]    -    m(x)E[en;T       <    S,  ] 

X        _L  X        J.  j.        X    —        -L 


=    -m(x)E[e1;Tx    <_  S±]  . 


Thus, 


NUM|2    =    {m(x)  (E[£l;Tx    <_    S±]     +    Efe^-S,     <    T    ])  (2.12) 

+    ?2E[Tx;Tx    <    S1]}2    +    ^2(-m(x)E[ei;Tx    <    Sj)2- 


From   Lemma    1    on    page    512    of    Feller    [1971] 


E[|El|;Tx  <  Sl]    <  |<5f4-,2l!lTx'Tx  isl]  (2-13) 


Et|c2|,Sl    <    Tx]     <    I(^|r)2E[S2;S1    <    Tx]      .  (2.14) 


Thus 


NUMI    ll<ET50ME|Tx;Tx    ISJ     +    E[S2;Sl    <    Tx])  (2.15) 

+    C2E[T     ;T       <    SJ     +     kl     -|-T   ^    E[T2;T       <    S,  ]  . 
xx—      1  '     '    m(x)    2  x      x  —      1 


Further, 


c-i^d|     >     |C|    E[S1;S1     <    T    ]     .  (2.16) 

—  xxx 


Thus 


a  +i^    -    .    -1.  JE1S,  ;S1     <    T    ] 
c  -  i£d        1  -  2.E,  '         11  x 


(2.17) 


±5^1    (E[Tx'*TxlSl]    +    E[S1;S1    <    Tx]) 

r2     i 

+     U|     E[Tx;Tx    <    SJ     +    ^-EtT^    <    S^  . 


Application   of   Schwartz'    Inequality    yields 


h  '^^ 


r=nl  E[si;Si  <  Tx] 


(2.18) 


<    o    1    v     E[S2]     +     (1    +    J  f  L  )E[T    ;T       <    S,  ] 
—   2m(x)  1J  2m(x)  x      x   —      1 


2m 


w  EISi]  +  (1  +  atfir'  <Eisii1/2  p(Tx  i  si)1/2» 


Let    3(x)    = 


P{T      <    S, } .       Then 
x   —      1 


m(x)    = 


E[min(Tx,S1) ] 


(2.19) 


Thus  , 


(rniW"  r^u1  E[si;Si     Tx] 


(2.20) 


<    3(x) 


E^S1^  2    1/2 

2E[min(Sn ,T    ) ]     +    E[Sl] 


6(x) 


1/2 


3(x) 


+  u 


3/2e!s2J1/2 


2E[min(S, ,T    ) ] 


THEOREM  1.   Under  Assumptions  (2.  5) -(2. 6) 


T 

p{mT$r  -  y}  "  (1  ~e~Y) '  -  0(p{Tx  -  si}1/1°) 


Proof.   It  follows  from  Lemma  1  that 


T 

^f1    elC|  ElexpliCj^-}]     -     (1-iC)"1)  |E[S1;S1    <    Tx]  (2.21! 

E  [S    ] 
1    l5H6<x»    2Etmin(Sl,T    )]    +   EIsl)1/2    6(x)1/2) 

1         X 

E[S2]1/2 
+     151        <B(X)  2E[n,in(Sl,T    )]' 

1         X 


where    3(x)    =    P{T      <_   S,  }    as    before. 

Thus,    applying    Lemma    2    on    page    538    of    Feller    [1971]    for 

T    =    3(x)~1/5, 


tt|p{T       <    y}    -     (1  -e    Y)  I  243  (x)17" 


7 
E[SZ 


+   E[Sl;Sl    <    Tx]{3(x)-2/5[6(x)    2E[mi;(s      T    )}    +    E[S2]1/23(x)V2 

_L    X 

E[S2]1/2 
+  B(x)-3/5  |  B(x)3/2        X 


3  M  v  '  2E[min(S,  ,T  )  ] 

_L    X 


The  result  now  follows 


EXAMPLE.   In  Example  1 


x  °° 

(x)  =  F(x)  +   /   F(dy)   /   H(dt)  [1  - exp{-y (y ) G (y , x-y ) t } ] 
0  0 


x 

=1-   /   F(dy)h(y(y)G(y,x-y) )  (2.22) 

0 


~  —  st 

where  h(s)  =   J   e    H(dt)  is  the  Laplace  transform  of  H  and 

0 

F(t)  =  1  -  F(t).   If  C  and  S  are  two  independent  random  variables 

having  distributions  F  and  H  respectively,  then 


$(x)  =  1  -  E[exp  {-  y  (C)G(C,x-C)S};  C  <_   x]  .  (2.23) 

Consider  two  independent  load  combination  processes 
{X1(t);t    0}  and  {X  (t);t  >  0}  of  the  type  in  Example  1. 
Assume  that  the  conditional  distribution  of  the  shock  load  process, 
given  the  constant  load  process,  is  the  same  for  both  load  combinati 
processes;  the  constant  load  magnitudes  have  distributions  F, 
and  F~  respectively  with  F,  (t)    F-(t);  and  the  times  between 
constant  load  changes  have  distributions  H,  and  H„  respectively 
with  H,  (t)  <_   H~(t)  .   Let  C.  and  S.  have  distributions  F.  and 
H.  respectively  for  i  =  1,2.   If  y  ->  y  (y )  G  (y  ,x-y)  is  an 
increasing  function  of  y  for  y  <_   x,  then 

-y (C  )G(C, ,x-C  )S 
61(x)  =    1  -  E[e     L.  1  L       1;C1  <  x] 


10 


-p(C  )G(C  ,x-C  )S 
3]_(x)  <  1  -  E[e     Z     ^     A       Z;C2    <_   x] 


=  32  (x)  .  (2.24) 


Thus,  Theorem  1  suggests  that  the  convergence  to  exponential 
of  the  distribution  of  the  time  for  the  load  combination  process 
to  exceed  a  level  x  is  faster  for  load  combination  process  1 
than  that  for  process  2.   This  behavior  has  been  seen  in  the 
simulation  studies  reported  in  Section  4. 


3.   THE  TAIL  OF  THE  DISTRIBUTION  OF  T   FOR  FINITE  x 

x 


In  this  section  we  will  study  the  behavior  of  P{T   >  t} 

2  x 

for  finite  x  and  large  t.   The  probability,  P{T   >  t} , 
satisfies  the  following  renewal  equation 


t 

P{T   >  t}  =  PlS,  >  t,T   >  t}  +   f   L  (du)P{T   >  t-u}   (3.1) 
x  1       x         0  X 


where 


L  (t)  =  P{S,    t,T   >  S,}  .  (3.2) 

x  1  —     x     1 


We  will  assume  L  (0)  =  0. 

Following  the  argument  on  page  376  of  Feller  [1971]  we 
will  assume  that  for  each  x  there  exists  a  constant  k(x)  such 
that 


11 


/   eK(x)u  L  (du)  =  1  ,  (3.3) 

2\ 


Y(x)  =   /   teK(x)t  L  (dt)  (3.4) 

0  x 


is  finite,  and  the  function 


g  (t)  =  eK(x)t  P{Sn  >  t,T   >  t}  (3.5) 


is  directly  Riemann  integrable. 

It  now  follows  from  the  Key  Renewal  Theorem  that 


lim  eK(x)t  P{T   >  t}  =  -^y   /   g  (t)dt  .  (3.6) 


Since  g  is  assumed  to  be  directly  Riemann  integrable, 
integrating  (3.3)  by  parts  yields 


1  =  P{T   >  S  }  +  <(x)   /   eK(x)t  P{Sn  >  t,T   >  S.}dt.  (3.7) 
X      _L  -  J.         x      X 


Thus 


P{T  <_   S  } 
k(x)  =  — — ^ ± .  (3.8) 


/   eK(x)t  p{s,  >  t,T   >  Sn  }dt 
'  lxl 


12 


It  follows  from  (2.4)  that 


E[min(T  ,S, ) ] 
K(x)m(x)  =  — ± .        (3.9) 

/   eK(x)t  P{S,  >  t,T   >  S,  }dt 
0  ixi 


The  defining  equation  for  k(x),  (3.3),  (3.8),  and  (2.6) 

imply  that  as  a  function  of  x,  k(x)  is  nonnegative  and  decreasing 

with  lim  k(x)  =0.   If  it  is  further  assumed  that  there 

x+°°  0S, 

exists  0  >  0  such  that  E[e    ]  <  °°,  then  it  follows  from  the 

dominated  convergence  theorem  and  (3.9)  that 


lim  k (x)m(x)  =  1  .  (3.10 

x^°° 


EXAMPLE.   In  Example  1,  assume  that  H(t)  =  e    , 
u(y)  =  y,  and  G(y,x)  =  G(x).   Then  k (x)  satisfies  the  equation 


1  =   /   F(dy) — * (3.11) 

0  X  +  yG(x-y)  -k  (x) 


where  G(x)  =  1  -  G(x) . 
Further, 


,.,     fX   <(x)t   -  t   -  G(x-y)t  „,«. 
g  (t)  =   J   e       e     e       J         F(dy 

X       0 


Thus 
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x  1  1 

/   g(t)dt  =   /   F(dy)  ^ =  y  •     (3.12) 

0  0  A  +  yG(x-y)  -  k  (x)    A 


Therefore, 


limeK(x)t  P{Tx  >  t)  -^  (3.13) 


where 


oo      x  — 

i    \     r    r   «/j  x     k(x),  -At   -yG(x-y)t  ,, 
y(x)     J    J   F(dy)  t  e  v  ' A e     e       1         dt 

0    0 


=   /   F(dy)  — i .  (3.14) 

0  (A  +  yG(x-y)  -  K  (x)) 


If  it  is  further  assumed  that  A  =  y  =  1 ,  then  it  follows 
from  (3.9)  and  (3.11)  that 


x  l 

/   F(dy)  -±- 


k  (x)m(x)  =  °- l+G(x-y) (3.15) 


/   F(dy) 

0  1  +  G(x-y)  -k  (x) 


/   F(dy)  -± 

0  1  +G(x-y) 


Thus  k(x)  <  —. — r  in  this  case. 
-  m(x) 
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-x 


In  the  further  special  case  in  which  F(x)  =  G(x)  =  e   ,  X  =  1 
and  y  =  1,  <(x)    satisfies  the  equation 


x 
1  =  f (k(x))  =   /   e" 
0 


■y 


l+e-^-yJ  -k(x) 


dy 


(3.15) 


where 


f(l)  =  ^[eX  -  e  X; 


3.16) 


-x 


and  for  y  ^  1  such  that  0 . <  y    1  -  e   , 


f(y)(l-y)2  =  (l-y)d-e  x)-e  X[ln((l-y)+e  x)-ln((l-y)e  X+e  x) ] 

(3.17) 


Further , 


Y(x)  =   / 


0    (1  +e  (x  y)  -k  (x))2 


e  ^  dy ;   or 


Y(x) (1-k (x) ) 3  =  (l-K(x))(l-e  X)-2e  X [ In ( ( 1-k  (x)  ) +e  X) 


-  ln( (1-k  (x) )e  x+e  X) ] 


(3.1 
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4.   A  SIMULATION  STUDY  OF  THE  ACCURACY  OF  THE  EXPONENTIAL 
APPROXIMATIONS  TO  THE  FIRST  PASSAGE  TIME  DISTRIBUTION 


In  this  section  some  results  of  a  simulation  study  of  the 
accuracy  of  using  the  asymptotic  results  of  Sections  2  and  3 
to  approximate  the  distribution  of  T   are  reported.   All 
simulations  were  carried  out  on  an  IBM  3033  computer  at  the 
Naval  Postgraduate  School  using  the  LLRANDOM  II  random  number 
generating  package  (see  Lewis  and  Uribe  (1981)). 

The  model  of  Example  1  was  simulated  for  various  cases 

of  distributions  H,  F,  G  and  shock  arrival  rate  \i .      Each 

realization  simulated  the  sample  path  of  the  process  (X(t); 

t    0}  and  the  first  passage  times  0  £  T    ^  T    <_   .  .  .  £  T 

12  n 

were  recorded  for  several  levels  x,  <  xn  <  . . .  <  x  .   The  number 

1  —   2  —     —   n 

of  replications  was  5,000.   Sample  moments  and  quantiles  were 

computed  for  each  T   .   A  more  detailed  account  of  the  simula- 

i 
tion  can  be  found  in  Noh  (19  84). 

4.1   The  Exponential  Approximation  P{T   >  t}  =  exp{-E[T  ] t} 

a.   Model  A 

Tables  1-1,  1-2  and  1-3  report  the  simulation  results 

for  the  model  in  which  H  =  F  =  G  are  all  unit  exponential  and 

\i         1.   Table  1-1  reports  simulated  sample  mean  and  coefficient 

of  variation  for  T   for  various  levels  of  x.   As  expected, 

x  ^ 

the  simulated  coefficient  of  variation  approaches  the  exponen- 
tial distribution's  value  of  1  as  x  gets  larger.   To  assess 
the  quality  of  the  exponential  approximation  1  -  exp{E[T  ]   t} 

to  the  distribution  of  T  ,  quantiles  from  the  simulated  data 

x   ^ 

were  computed  and  compared  to  the  approximating  exponential 

16 


TABLE  1-1 
Simulated  Moments  for  T   in  the  Case  F  =  G  =  H  =  exp(l)  and  y  =  1 


x 


x-Level  E [T  ]  Coeff  Var [T  ] 


0.50  0.276  1.984 

M0.01) 


1.00  0.669  1.451 

(0.014) 


2.00  2.000  1.177 

(0.033) 


3.00  4.878  1.080 

(0.074) 


4.00  11.539  1.024 

(0.167) 


5.00  26.945  1.002 

(0.382) 


* 


(  )  is  the  standard  error  of  the  mean 
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quantile  q 

a 


q   =  -E[T  ]  ln(l-a)  (4.1) 

a       x 


where  for  this  model  E[T  ]  can  be  computed  analytically  and  is 
found  in  Gaver  and  Jacobs  [1981].   Table  1-2  reports  the 
simulated  quantiles  with  the  approximating  exponential  quantiles 
below  in  parentheses.   As  expected,  the  exponential  approxi- 
mation is  better  for  the  large  level  x  =  5  than  for  x  =  0.5. 
The  approximation  is  also  better  for  a  >  0.5. 

One  way  the  distribution  of  T   differs  from  an  exponential 
is  that  it  has  an  atom  at  0;  in  particular  for  the  model  A 
p{T   =  0}  =  1  -  F(x)  =  e   .   The  sample  quantiles  for  the 
simulated  conditional  distribution  of  T   given   T   >  0  appear 
in  Table  1-3.   Below  the  simulated  quantiles  in  parentheses 
appear  approximating  exponential  ones  computed  as 

-E[T  ] 
%    '    PIT   >  0)  ln(1"a)  •  (4'2) 

X 

The  quantiles  of  the  simulated  conditional  distribution 

are  much  closer  to  their  exponential  approximation  than  those 

for  the  unconditional  distribution, 
b.   Model  B 

Tables  2-1,  2-2,  and  2-3  report  simulation  results  for 
the  model  which  is  the  same  as  Model  A  except  that  the  distri- 
bution of  the  constant  load  magnitudes  F  is  exponential  with 
mean  1/2.   The  simulated  means  and  coefficient  of  variations 
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TABLE  2-1 
Simulated  Moments  for  Model  B 


Level  x  E [T  ]  Coeff  Var [T  ] 

x  x 


0.5  0.577  1.453 

*(0.010) 


1.0  1.329  1.144 

(0.022) 


2.0  3.980  1.036 

(0.058) 


3.0  10.530  1.006 

(0.150) 


4.0  28.014  0.997 

(0.395) 


5.0  75.425  1.003 

(1.070) 


(  )  is  the  standard  error  of  the  sample  mean 
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appear  in  Table  2-1.   An  analytical  expression  for  E [T  ]  for 
this  model  appears  in  Gaver  and  Jacobs  [1981] .   This  expression 
is  used  in  the  approximating  quantiles  (4.1)  and  (4.2).   The 
approximating  quantiles  for  the  conditional  distribution  of  T 
given  T   >  0  which  appear  in  Table  2-3  are  closer  to  their 
corresponding  simulated  quantiles  than  the  approximating  quan- 
tiles for  the  unconditional  distribution.   Comparison  of 
Tables  2-2  and  2-1  for  the  x  >_  2  and  a  <_   0 . 3  suggests  that  the 
distribution  of  T   is  converging  faster  to  an  exponential  for 

Model  B  than  for  Model  A.   It  follows  from  (2.2)  that  P{T   <  Sn } 

x     1 

is  smaller  for  Model  B  than  for  Model  A.   Thus,  Theorem  1 

suggests  that  the  convergence  of  the  distribution  of  T   to 

exponential  should  be  faster  for  Model  B  than  for  Model  A. 

1    -k (x) t 
4.2   The  Exponential  Approximation  — ; — r-  e 

f  rr  Y (X) 

In  this  subsection  simulation  will  be  used  to  study  the 
exponential  approximation  suggested  by  the  asymptotic  result 
(3.6).   This  is  an  approximation  for  P{T   >  t}  for  fixed  finite 
x;  it  should  be  more  accurate  for  t  large. 

Two  cases  of  Model  1  were  simulated.   In  both  cases,  shock 
loads  arrive  according  to  a  Poisson  process  with  rate  1  and 
constant  loads  change  magnitude  at  the  times  of  arrival  of  a 
Poisson  process  with  rate  1;  the  shock  level  magnitudes  have 
an  exponential  distribution  with  mean  1.   In  Case  A,  the  dis- 
tribution of  the  constant  load  magnitude  is  exponential  with 
mean  1;  in  Case  B,  it  is  exponential  with  mean  1/2. 

In  both  cases  considered,  it  is  possible  to  determine 
analytical  expressions  for  the  integrals  determining  k(x)  and 
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y(x).   The  value  for  k(x)  was  found  numerically.   Values  of 
k(x)  and  y (x)  for  the  two  models  can  be  found  in  Tables  (3-1) 
and  (4-1).   Note  that  as  x  increases  k(x)  decreases  and 
approaches  E[T  ]   .   As  x  increases  y(x)    decreases  and  approaches 
1.   As  expected  k(x)  £  E[T  ]~   for  all  levels  of  x. 

To  assess  the  accuracy  of  the  exponential  approximation, 
quantiles  of  the  simulated  data  were  computed .   These  quantiles 
appear  in  Tables  (3-2)  and  (4-2).   For  each  level  x,  the  first 
row  gives  the  simulated  quantile,  the  second  row  gives  the 
approximating  exponential  quantile 


Q*  =  "  ^y  ln(Y(x)  (1-a))  ;  (4.3) 


and  the  third  row  gives  the  approximating  exponential  quantile 


Q+  =  -E[T  ]  ln(l-a)  .  (4.4) 


The  exponential  approximation  (4.3)  is  in  general  closer 
to  the  simulated  quantile  than  (4.4).   However  the  two  approxi- 
mations become  closer  as  x  gets  larger.   As  expected  (4.3) 
approximates  well  the  simulated  quantile  Q   for  a  >  0.75  for 
all  values  of  x.   However  if  x  is  sufficiently  large  (4.3) 
approximates  fairly  well  the  simulated  Q   for  a  as  small  as  0.1. 
A  comparison  of  Tables  (3-2)  and  (4-2)  suggests  once  again  that 

the  convergence  of  the  distribution  of  T   to  exponential  is 

x 

faster  for  Model  B  than  for  Model  A. 
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TABLE  3-1 
Values  of  k(x)  and  y {x)    for  Model  A 


Level    x 

K(X) 

ElTx] 

Y  (x) 

0.2 

1.708 

9.492 

5.942 

0.4 

1.459 

4.486 

3.254 

0.6 

1.247 

2.814 

2.364 

0.8 

1.066 

1.977 

1.923 

1.0 

0.912 

1.475 

1.670 

1.2 

0.780 

1.142 

1.502 

1.4 

0.667- 

0.906 

1.386 

1.6 

0.570 

0.730 

1.303 

1.8 

0.487 

0.596 

1.238 

2.0 

0.417 

0.491 

1.193 

2.2 

0.355 

0.407 

1.155 

2.4 

0.303 

0.339 

1.125 

2.6 

0.259 

0.284 

1.102 

2.8 

0.221 

0.238 

1.085 

3.0 

0.188 

0.201 

1.068 
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TABLE  4-1 
Values  of  k(x)  and  y(x)  for  Model  B 


Level    x 

k(x) 

E[Txrl 

Y  (x) 

0.2 

1.563 

4.760 

3.109 

0.4 

1.237 

2.265 

1.860 

0.6 

0.975 

1.433 

1.400 

0.8 

0.801 

1.015 

1.286 

1.0 

0.650 

0.762 

1.183 

1.2 

0.531 

0.592 

1.122 

1.4 

0.435 

0.470 

1.084 

1.6 

0.358 

0.378 

1.059 

1.8 

0.296 

0.306 

1.044 

2.0 

0.244 

0.250 

1.032 

2.2 

0.199 

0.204 

1.019 

2.4 

0.166 

0.168 

1.017 

2.6 

0.136 

0.138 

1.011 

2.8 

0.112 

0.113 

1.007 

3.0 

0.092 

0.093 

1.005 
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